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Abstract 

It has been shown by Son and Surowka that the presence of anomaly in hydrodynam- 
ics with global U{1) symmetry can induce vortical and magnetic currents. The induced 
current is uniquely determined by anomaly from the existence of an entropy current 
with non-negative divergence. In this work, we extended the analysis to hydrodynamics 
with [/(I) symmetry spontaneously broken, i.e. U{1) superfluid hydrodynamics. We 
found that all possible first order gradient corrections are determined up to five arbitrary 
functions, with the entropy current containing one arbitrary function. Furthermore, the 
stress tensor does not receive correction from terms proportional to the magnetic field. 
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1 Introduction 



The hydrodynamics is beheved to be a universal description of quantum field theory in 
long time and large distance limit. The assumption of local thermal equilibrium allows 
a description of the system in terms of fluid velocity u^{x) and local thermodynamical 
quantities including temperature T{x) and chemical potential fi{x) for conserved quantities, 
which has the slowest relaxation. If the underlying quantum system has a spontaneously 
broken symmetry, the resulting gapless Goldstone mode should also be taken into account 
as an additional degree of freedom. Normal hydrodynamics should be replaced by superfluid 
hydrodynamics, the dynamics of which includes a superfluid component with velocity ^^(x). 
There have been significant efforts in formulating superfluid hydrodynamical description of 
QCD with spontaneous chiral symmetry breaking[l, 2, 3]. Superfluidity is also believed to 
be relevant for the phenomenology of heavy ion collisions [4]. 

When the system under consideration is not evolving slowly or smoothly enough, 
the ideal hydrodynamical equations should be improved by a systematic inclusion of of 
gradient terms, compatible with the symmetry of the system. The coefficients appearing 
in front of the gradient terms are referred to as transport coefficients. For normal hydro- 
dynamics, the transport coefficients corresponding to first order gradient terms are shear 
viscosity, bulk viscosity and charge diffusion constant. They should be calculated from the 
microscopic theory. For example, in deconfined phase of QCD or QCD-like theories, the 
transport coefficients have been extensively studied in both weak coupling[5] and strong 
coupling regime[6, 7, 8, 9]. For superfluid hydrodynamics, the number of transport coeffi- 
cients increases significantly due to the additional superfluid velocity. It has been recently 
enumerated in [10], see also [11] that there are 14 transport coefficients for U{1) superfluid. 
The transport coefficients have also been calculated by method of gauge/gravity duality. 

While previous studies have been focusing on the parity even term in the gradient 
expansion, parity odd terms can be equally important. Early studies of fluid/gravity duality 
found an unexpected term proportional to the vorticity of the fluid in the constitutive 
equation of the R-current[12, 13]. The puzzling term seems to lead to negative entropy 
production at the first sight. This issue is clarified by Son and Surowka in [14], where they 
showed that the vortical term, as well as a term proportional to the magnetic field is actually 
required by the existence of entropy current with non-negative divergence. Working to the 
first order in gradient expansion, they found that in the presence of global U{1) anomaly, 
both the charge current and the entropy current get modified by a vortical and a magnetic 
term. Remarkably, the form of the correction is almost entirely determined by the anomaly. 
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This result was generalized to nonabelian symmetry group in [15], where similar effect has 
been found. 

In this work, we would like to consider the effect of triangle anomaly on the U{1) 
supcrfluid hydrodynamics. Wc will focus on the parity odd sector and use the existence of 
entropy current to determine the transport coefficients. We will first write down the ideal 
superfiuid hydrodynamical equations in external U{1) field in Section 2. We will construct 
all possible first gradient terms in parity odd sector in Section 3, following the technique 
in [10]. The constraint equations on the transport coefficients will be derived based on the 
criterion of non-negative entropy production in Section 4. The constraint equations are 
solved in Section 5 and we will conclude in Section 6. 

2 Ideal superfiuid liydrodynamics in the presence of external 



The nonrelativistic hydrodynamical equation for C/(l) superfiuid in the absence of external 
field is known as Landau's two fluid model[16]. It can be derived with different approaches, 
including quasi-average[17], Poisson bracket[18] and effective Lagrangian[19, 20]. The gen- 
eralization to relativistic superfluid hydrodynamics is written down in [21, 22, 23, 24, 25]. 
We quote the result of ideal relativistic superfluid hydrodynamics of [25] in the following: 



The first and second equation in (1) is energy-momentum conservation and current 
conservation respectively. The third equation is a "Josephson" equation. It is a statement 
that the density of the fiuid and phase of the condensate is not independent, with (p being 
the phase of the condensate. The constitutive equations are: 



field 



< 



d^j^ = 

^ /Lt + u'^dij.ip = 0. 



(1) 



(2) 
(3) 



together with the equation of state: tie = Tds + /idn + f ^dX and dp = sdT + 
nd^ — f^dX. Here X = ^{dip)'^ is the additional thermodynamical variable due to the 
superfluidity. The entropy is conserved in the ideal superfiuid hydrodynamics: 
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d^isu^") = 0. 



(4) 



In the presence of an external U (1) field (provided that it does not destroy the superflu- 
idity), the hydro dynamical equations need to be modified. As in normal fluid, external elec- 
tromagnetic field will exert Coulomb and Lorentz force on the fluid. The energy-momentum 
is not conserved as in the isolate system of superfluid. A term of F^^j\ is needed on the 
RHS of (1) to account for it. Furthermore, assuming the external field couples to the 
matter field through minimal coupling in the microscopic theory, the superfluid velocity re- 
ceives contribution from the external field in the following way: — >■ V<p — A. Taking into 
account relativistic invar iance, it is not difficult to convince oneself that the hydrodynamical 
equations of superfluid are given by: 

daf = 

< , (5) 

with Tf"" = {e + p)ui^u'' + pr]!^" + P^^C" and j" = nu" + p^f". The equation 
of state is the same as the case without external field: de = Tds + ndn + f^dX and 
dp = sdT + ndjj, — f^dX, but with X = The phase of the condensate if is replaced 
by ^ju, the un-normalized superfluid velocity. We have confirmed the intuitive equation (5) 
with Poisson bracket method [31]. Starting from ideal hydrodynamical equations, we can 
show that there is a conserved entropy current: The first two equations of (5) give, 

^u.d,T^^ + !^d,f = ^Ur^F^^h. (6) 
Upon using 5,/^^ = —F^j^i, and the equation of state, we readily obtain dn{su'^) = 

0. 

3 First order gradient expansion in the parity odd sector 

In this section, we want to study the effect of axial anomaly on superfluid hydrodynamics, 
which amounts to a modification to the current conservation equation: d^jf^ = CE B, with 
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_ and = ^e^^'^^UnFa^. In this work, we will restrict ourselves to a single 

L'"(l) charge. The anomaly can be realized as an AAA anomaly in a fictitious U (1) theory: 
C = ■07'^((?^ + ie7^A^)^, with -ip and being the matter field and external axial U{1) field 
respectively. It will be interesting to investigate the case of two charges for phenomelogical 
applications [27, 28, 29, 30]. We leave it for future work. Note that the superfluidity arises 
from spontaneous symmetry breaking. We should ensure that the ^7(1) symmetry (charge 
conservation) is not spoiled by the anomaly term CE • We can for example turn 

off the electric field, keeping only the magnetic field, so that the anomaly term vanishes. 
However its effect on charge current and entropy current remains, as is clear from [14]. 

The first order gradient expansion will introduce correction to the stress tensor, cur- 
rent and chemical potential[ll, 10]: 



T"" = (e + p)u^'u'' + pr)>"' + f^^" + tt"" (7) 
f = nu^" + f^^" + u^" (8) 
/i = ■ u + /XA- (9) 

We can specify a frame to further constrain the corrections. The frames often used 
are the fluid frame, where tt^uU^ = and = and the transverse frame, where tt^uU^ = 
and I'fxU^ = 0[10]. We will choose to work in the transverse frame in the following. 

The next step is to to write down all possible parity odd, first order gradient terms 
that can appear in tt^j,, i/^ and /xa- These include gradient of the fluid velocity (Ou^ and 9^^) 
and gradient of the thermodynamical variables and dX). However it is complicated 

by the fact that they are not completely independent, but related by ideal hydrodynamical 
equations. The procedure we will use in enumerating independent gradient terms closely 
follow the method adopted in [10]: We flrst list all possible (pseudo)scalars, (pseudo)vectors 
and (pseudo) tensors, which are independent upon using the ideal hydrodynamical equations. 
Then we construct tt'^'^, and ha out of them. 

Assuming the external U{1) fleld ~ 0(p^), the fleld strength = dij,A,y—di,An is 
first order in gradient. There are in total 2 pseoduscalars, 7 pseudovectors, 2 pseudotensors, 
7 scalars, 7 vectors and 2 tensors that are first order in gradient. Note that due to the 
external field the number of independent structures are slightly larger than those reported 
in [10], where the numbers are 1, 5, 6, 5 for pseudoscalar, pseudovector, scalar and vector 

*^Thcrc is also an interesting effect due to gravitational anomaly, see [26]. We do not discuss it here 
*^We thank Ingo Kirsch for pointing this out to us 
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respectively. The independent structures are listed as follows, 
pseudotensor: 



X"" = e^'^'^Upi^Xl + e'<"'^Upio.\% (10) 



pseudovector: 



= l^e^'^'^u^iadpu^, (11) 
pseudoscalar: 

i-^.i-B, (12) 



tensor: 



U 



vector: 



X^.v - -P^ I 2 Va^—^ j (I'Jj 



(T^,J'', d^p, d^n, d^X, i,d^u\ E^, F^^r (14) 



scalar: 



l^d^^p, l^'d^X, u^d^^X, E ■ (15) 

where P^'^ = ij^'^ + u^u'^ is the projection operator, which is transverse to the normal fluid 
velocity u^P^'^ = 0. l'^ = P^^^,i, is the component of transverse to the normal fluid 
velocity. 

Now we are ready to decompose ir^^ , u^^ and in terms of the independent struc- 
tures. Note the constraint u^'k'^'' = 0, u^u'^ = from choosing the transverse frame, we 
have: 
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+ UkI^^'K''^ + Uul^^W"^ + Uvl^^V^ + HlZ^'^L^)) + l^'riN^^ -u + Nb^-B) (16) 

1^" = {V^co" + Vb^^ + Vn^f" + y^iC'^ + VuU" + "IV^'' + Vl^^'') + (A^C • w + Afi^ 5)^'' 

(17) 

liA = A^^-oj + Ab^ ■ B, (18) 
where A'-i^B"^ = \ {Ai^B" + A'^B'^). 

4 Constraint from entropy conservation 

In the presence of triangle anomaly, the charge current is no longer conserved: d^j^ = 
CE B. This enters as a second order correction in gradient to the hydrodynamical equations, 
thus wc need to use the constitutive equations up to first order in gradient expansion, which 
has been worked out in the last section. Therefore, we have: 

daf = CEB 
< , (19) 

lJi + HA = -i-u 

Note an explicit change to the form of "Joscphson" equation due to the gradient correction. 
This is because wc will choose as a separate thcrmodynamical variable, which is always 
zeroth order in gradient. While the gradient corrections to the energy-momentum and 
current conservation are implicit through the constitutive equations: T^'^ = (e + p)u^u^ + 
Pf]^^ + /^C^C^ + T^^" and = nu^ + /^^^ + u^^. In the presence of first order gradient 
correction and anomaly, the entropy current defined in (6) is not conserved anymore. It 
satisfies the following equation: 

d,{su^ - f^u^) = -^d.u.TT'^'' - d, (I) u'^ + ^d,{nu^) + ^-t.CE.B. (20) 

We should modify the definition of the entropy current and the positivity of entropy 
production will allow us to fix the transport coefficients. The parity even sector has been 
investigated in [10, 11]. In particular, this criterion can significantly reduce the number of 
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nonvanishing transport coefficients. We are interested in the parity odd sector. Remarkably, 
in case of normal hydrodynamics, Son and Surowka showed the positivity (rather vanishing 
for the parity odd sector) allows them to almost completely fix the transport coefficients. 
As we will sec soon, this is also true for supcrfluid hydrodynamics. 

We start by writing down the most general form of the entropy current, up to first 
order in gradient: 

s/* = su^" - ^1^" + s^u'' + sbB" + snQ" + skK" + suU" + syV + slLI" 

+ {Snuj^ ■ + SnB^ ' B)l^' + {s^^ ■ CO + SuB^ ' B)u^' . (21) 

The corresponding entropy production is given by: 

sn^" + skKI" + suU^" + svV + slL" + {s^^^ ■ u + • S)/'^ + (s„^^ • oj + s^b^ ■ B)u^'). 

(22) 

We first note that all terms in the expansion of (22) are parity odd, second order in 
gradient. Most of them can be decomposed on a basis in the form: pscudoscalar-scalar, 
pseudovector-vector, and pseudotensor -tensor, with the independent structures listed in the 
previous section. However, the divergences of the last four terms of the entropy current con- 
tain new structures with double derivatives: l^^^d/^iUj^, l^^^^d^B^,, u'^^'^d^ujy and u'^^''d^B„. 
Only three of the four are independent of each other. The last two are related by ideal 
hydrodynamical equations. We start with the following identity from ideal hydrodynamical 
equations: 

+ p)u^'u,) + d,p + fd^X - C,5^(nw^) = nE,, (23) 

where w = e + p is the enthalpy. Applying e^^'^^^aUpdx to the above, we obtain: 

2^ • u^|,{wu^') - 2Q^'df,w + 2w{i ■ uO - Lo^l^d^u^) + we^"''''^^^u^,uPdpdaUfs + C ■ Bd^inui") 
- e^'^'^^d^fu.^^dpX - Ki'd^n + ne'^'^^^i^u.uPdpdaAp - niBn'^d^u, - ^ ■ BO) = 0. 

(24) 

In addition, snd^^^^^ and SKd^K^ also contain the following contributions to the structure 

ui'i^d^.uj^ and uf'^^'d^B^: 
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The requirement of the vanishing entropy production on the three independent basis forces: 
Snw = SnB = 0, while the other two entropy current coefficients are related by: 

Suui = 'iwsu + sn SuB = nsu + sk- (25) 

We can then proceed to decompose all the terms on the basis given by pseudotensor 
•tensor, pseudoscalar-scalar and pseudovector-vector. In achieving this goal, we will need 
the following identities: 

a^u^n"- = a^„i,"- + ^4 (26) 



jf- (28) 



dp djjL dX 

It is easy to obtain from dp = sdT + nd^ — f^dX that 

d{ii/T)_^_j^ dj^/T) _ w wri_ 

dp sT^ dfi sT^ dX sT^' 

The last term d^{soOi'){0 = co, B,n, K,U,V, L) generates Oi^d,^{so) and sod/^O". The 
former can be brought into the form vector-pseudovector via 

Mso) = ^a,P+'-^a,, + ^a,x. (29) 

while the latter needs some extra work. By ideal hydrodynamical equations, d^O^ can 
be expressed as a linear combination of pseudoscalar-scalar, pseudovector-vector and pseu- 
dotensor -tensor. The only exceptions are d^^l^ and d^K^, which contain terms with double 
derivatives and arc used to fix the entropy current coefficients above. We listed the final 
results below and move the details to the appendix. 
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^^UJ^' = l.(nE-co- cu^d^p - fuji'd^X + ^ ■ ojQ (30) 

dr.Bi' = -(nE-B- Bf'd^p - fBi'daX + BO -2E -oj (31) 

w 

d^n" = -ei-uj- ^{-B^d^p - fB^d^X + nE-B + ^-BQ- ^{nE ■ co + ^ ■ uC 

- u^d^p - fco^d^X) + u'^^'^d^u, ~ ^e^^'^^C^.u^uPdpd^Uf} (32) 
d^K^" = -ei-B + 2iiE-u + ^{-Kf'd^p - fK^^d^X) + ^E ■ B - B'^d^./j. 
-e'^'^-^C^u^ufdpdaAf, (33) 
d^.U^' = ^{nKi'd^p - 2fU^'^^,X) - ^iu^'^^p - ^ ■ cou^d^p - ^B^'^^p (34) 
^^^V^ = ^ {nK^'di.fi + 2U^'d^fi - 2fV^d^X) - fiuj^'d^n - ^ ■ uui'd^fi - ^^Bi'd^ii (35) 

w 2 ft- w 

- fl'^d^X + fC) - \B>^£rd^u, + \d^ee^-"^u,edpux. (36) 



Note in the above, u^d^p, u^d^fi are not listed as independent scalars in the previous 
section. They can be expressed as: 

u^'^^f^ = bsi-Os) + 6„(C - nO) + bxuf'd^X (37) 

ui'di.p = Csi-es) + CniC- nO) + cxuf'd^X, (38) 



where 6. = (f) K = (g^)^^^ , 6. = {§^) and c. = (i) = ( 



s,X' 



\ / s,n 

Let us consider the basis of the form pseudotensor-tensor first. The only basis get 
popularized are cr^ya^^^ and a^j,y\^^^. The former receives the contribution from Tcra^^a^^^ , 
which vanished identically thus does not constrain Tcr[32]. On the other hand, the latter 
receives contributions from Txa^^uX'^'^ and slO^L'^. The last term of d^L^^ as shown in (36) 
should be expressed as follows: 



^crfj.u>^'^ +^-^{ ^^—^2 \-u'^dfj_ii+ 



nE.^-l^^d^p-pl^^d^X + l\ ^ 1 1 (f, ^<:-^ ^^ 
) - -i-BQ + -zOixvV'B^ + u^id^X + fx{df,fi + ^ dfj,Uu)) 

IV o Zi 

(39) 
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The vanishing of the component of entropy production on the basis a^i,\^^'^ gives: 

Tx = -jSL. (40) 

Now wc are ready to move to the basis of the form pseudoscalar-scalar and pscudovec- 
tor -vector. There is a last subtlety that not all pseudovector-vector are linearly independent. 
Some of them can be expressed as a linear combination of others and possibly also pseu- 
doscalar -scalar . We work out the relevant relations in the appendix. The complete basis 
set we will contains 2x7 pseudoscalar- scalar and 21 pseudovector-vector. The latter is 
listed as follows: 

c7^,/^5^ ^^r^^w., B^'Cd^u,, K^^Cd^.u,, uH^d^p, BH'^d^p, Kn^d^p, m^d^p, 

un^'d^ix, Bn^'d^ix, K^'rd^n, m^d/.n, u'^rd/.n, w'^ra^x, s^^ra^x, K^^ra^x, m^'d^x, 

Ui^rd^X, Vrd^X, E-uj,E B. (41) 

Plugging (36) to (22) and using (64) -(90) to organize all terms into the above basis, 
we obtain a set of 35 coupled equations upon demanding the non-negativity of the entropy 
production(Clearly all coefficients in the expansion on the basis have to vanish). Wc will 
not spell out the full equations, but only comment on one property of the equations. They 
are all linear in the unknowns, algebraic for the transport coefficients and first order PDE 
for the entropy coefficients so{0 = co, B,Cl, K,U,V, L,u). This property will be crucial in 
solving the equations. 

5 Solving the coupled equations 

As is noted in the previous section that the coupled equations are first order PDE for the 
entropy coefficients, but are only algebraic for the transport coefficients. We can simplify 
the problem by solving the transport coefficients in terms of the entropy coefficients. This 
step is relatively straightforward. We end up with the following results: 
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SL = 0, Ub = o, nu = o 

f sT 
-Hi = n^, /xFl = — nn 
2 w 



liu = /illv 



2MSir - 2sB + 2srs„ + — ^ + (1 ^)sn 



w 



w 



sT"^ 2Tw ~ sT 

Vb = -T [ -— + (2 - —)SK + + Su - — 

\ T 2w WW 2w 

T{nbn + sbs) , nsu , fiVx Vu Vq n 

wihcr, - b^c,) ^" + ^ + ^ + 2T " 2;T^ + ^^^^ " 2^^^^" = ° 



-2TCn + 2^Cs 
w{bsCn - bnCs) 

^Vu ( 2m(1 + f&)Mn + sb,) - ^{nbn + sb,) 



sv wVu ijVv Vu . 



sT2 



+ 



Su + ( yft^X - ;^ ) sc/ 



Ws{bsCn - bnCs 

+ (—dp \odx] sn = 

\w wj-^ J 

Pb = -nT (1 + {ncn + scs)dp + {nbn + sbs)d^) s« - T ((nc^ + scs)dp + (n6n + sbs)d^) 
wiJ.{nbn + sbs) 



sT 

Nb = ^Ab 



-Ab 



Ab = - 
Ab = -T 



,T2 



w{l + iJ,bx) 

S/isk SB nfx 



{{dx + cxdp + bxdfj,)sK + n{dx + cxdp + bxdfj)su) 

-'^n , At 2 



r, H Su - 7^ + -7^(1"^ + wbn)AB + {Cndp + 6„9^) SK 

2w w w 2w sl^ 
+ n {Cndp + bnd^j) Su 



2w 



cxsu - bxsv 
IIl w 



sT 



w{l + nbx) 



sT 
+ —Uv 
2w 



2w {dx + cxdp + bxdf^) Su + {dx + cxdp + bxdf,) sq 



= -2wTsu + T{ncn + scs)su + T{nbn + sbs)sv - ^{nbn + sbg) ( A^^ + ^Ily 



sT 



sT, 



2w 



2wT {{nCn + scs)dp + {nbn + sbs)d^) Su-T ((nc„ + scs)dp + (re^n + sbs)dfj,) sq 



-r 



2Sa; // 



- 2/Lts„ - c„S[/ - H ^-^{l + 'wbn){A^ + —Uv) 



w 2 

+ 2w{Cndp + bnd^)Su + {Cndp + bnd^)sQ, 



2w 
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We find the transport coefficients are not fully determined by (42) in terms of entropy 
coefficients sk,sb,Si^,su,sv,Su- In particular, we can choose Ily, 11^^, and Vh as the free 
inputs, in addition to T^r found in the previous section. 

The entropy coefficients, independent of (42), satisfy the following equation: 



(43b) 



w 



sa - i^sv = { - - TTFf'^'^ -dn]Sij (43c) 



SQ. - jJiSv = -IJ.sk + 2 ^1 - - iJid^^ SB 



(43d) 



S€i= - ) + 2/^*^ ~ 2^-5 (43e) 



2C u? ( IsT 

sn = Tf- + ( 2 - -j^dx ] SB + ^iSK (43f) 



- {j2^^ - ^p) + ^ (72 - ^p) + ^ " ^^^) " ° ^^^^^ 

-2rsc„ + 2//SC5 - 2/L((nCn + scs) + 2u;(n5„ + s6s) /„ 1 w \ 

+ 5^ h — jo^x su 



Ws{bsCn - bnCs) \ 1^ l^f 

+ ( - 9p) 5y - [Jj^dx - i^p) = (43h) 

wT{nbn + sbs) + fiT{nCn + scg) + sT{—Tcn + ncg) (n sT 



(ft sT \ 
2 ~^ 'Y^^ ) ^^^^ '"^'"^ 



w{bsCn - hnCs) 

There arc seven unknown functions to be solved for from the nine equations. It is 
convenient to eliminate s„ from (43h) and (43i) to obtain: 

-2 (^^dx + lid^ SK + {j^dx + s^-{l- '^dx^ sy = 0. (44) 

This is an over-determined problem. In order to proceed, we first take a moment to prove a 
commutation relation between the following differential operators jzdx—dp, jidx+lJ-d^ and 
1 — ^dx- Bear in mind that we have chosen p, n, X as the independent thermodynamical 
variables, w, s, T, p contain implicit dependences on them through the equation of state. 
The dependences are subject to the Maxwell relations. Starting from dT = <^P~'^<^i^+f '^^ ^ 



s 

we have 
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Using the above relations, it is a short exercise to show the operators jzdx — dp, jidx + M^/* 
and 1 — ^dx commute with each other. 

To solve the equations (43a)-(44), we apply jidx — dp to (43c)-(43d) to obtain: 



^dx - d^ - ^dx^ - 2 (^j^dx - dp^ (^2 - ^dx^ sb + ^ 



Using (43a), (43b) and noting the commutation relation of the operators, we obtain 



j^dx - dpj SK = 0. (47) 

Note that U = /^^j we conclude that the solution is of the form sk = sk{T, ii). Plugging 
(43d) and (43f) into (44), we obtain: 

(l + M5^ + ^ax)sK = 0. (48) 

With the specific form sk = sk{T, h), we can solve the above by: 

« = ^. (49) 
A* 

The singular behavior of s;^ as ^ forces sk to be zero. To proceed further, we note 
from (43a) and (43b) that —^ + 2sb is also annihilated by jidx — dp. Therefore, it can 
be parametrized as: 

-'-^ + 2sB = -ss{T,^,). (50) 

Next, we plug (50) to (43c)- (43d) and (43e)-(43f), the resultant PDEs give the following 
unique solution: 

It can be verified that the rest of the equations are identically satisfied. We list our results 
on the coefficients of entropy current in the following: 

SK = 

2Cm^ 

Su, = 2llSB — 



3T 

1 

2Cii (n 



su = -2[ -p^dx - dp] SB 



sv = -^ + '2[-^dx + d^]sB 



sn = — + 21- -TT^dx SB 



T ' \ p 
Su = 0. (52) 
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Plugging (52) back to (42), we obtain the following nonvanishing non-entropy coefficients: 



Pb = Nb = Ab 

n,. 



-Hi 
2 



w 



Ab = - 



w 



+ — Hy 

2w 



dxSB 



2sT^ 



+ 



w{l + iJ,bx) 
iP + cx-^bx),^'^ 



P + CX- nbx 
sT 



{dx + cxdp + bxd^] 



Y^xsbJ 



Ul _ w 



(2 2n \ 

P<^ = i {ncn + scs){-- + 2Tdj,) + {nbn + sbs){— + 2Tdf,) \ 



sT 
1 



2 dxSB ] - -^i^^n + Sbs) 



+ —Uv 
2w 



(nCn + SCs) — — + [nbn + sbs)- 



sT 



sT 



A. 



2{nbn - Cn) ^ijlT 



w 



+ 2T{Cndp + bnd^)^ i^^dxSB 



w 2Cfi^ 

^ 1- [ — bn - Cn) 

6W II 



Ts 



{t' + Wbn) ( A^ + ^Hy 



- — 



2wn 

sT^ fCfj 



2w 

2sT sT. 



-2T{1- —)-^dxSB + 



w \ 




2nT 




w 





- -pOxSB 



dx - dr 



sT 
7 



dxSB 



w 



liVv - Vh = 2nT 



+ 

1 
w 



w 
1 

7 



sT 



sT 
P 



dxSB 



;dx 



d, 



w 



sT . 

7 



rdxSB 



(53) 



Summarizing the results, we have determined, within the parity odd sector, the first order 
gradient corrections to the stress tensor, charge current and "Josephson" equation up to 
five arbitray functions. The five functions are chosen as sg, Ily, Iltj, Vh and Tg-. In 
particular, the entropy current, as well as the "Josephson" equation is parametrized by 
the function sb, and the rest enters the stress tensor and charge current. It is interesting 
to note that terms proportional to the magnetic field is absent in the correction to stress 
tensor, i.e. Pb = Nb = Ab = 0- The presence of arbitrary functions in first order 
corrections is in contrast to the case of normal hydrodynamics [14, 15], where the positivity 
of entropy production almost fully constrain the first order gradient corrections. We can 
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also compare our results with [32] (see also a more recent generalization in [33]). The authors 
of [32] studied the constraint of the positivity of entropy production in first order superfluid 
hydrodynamics, with both parity even and odd terms. A total of 6 arbitrary functions are 

needed to parametrize the full gradient correction. 2 of them appear in the entropy current. 
The authors claimed a nontrivial mixing between the parity odd and even sectors. While this 
is true in general, the nature of anomaly that is temperature and density independent (see 
for example [34]), seems to suggest its induced effect to be non-dissipative(care is needed in 
a precise definition of anomaly induced terms). More recently the vanishing of the entropy 
production in the parity odd sector in the context of normal hydrodynamics has been argued 
in [35] based on the principle of time reversal invariance. It will be interesting to see whether 
it is also true in the case of superfluid hydrodynamics. 

6 Conclusion 

We have extended the study of anomaly effect on the constitutive equations of normal 
hydrodynamics[14, 15] to superfluid hydrodynamics. We have enumerated all possible 24 
parity odd, first order gradient corrections to the constitutive equations: 13 for the stress 
tensor, 9 for the charge current and 2 for the "Josephson" equation. We used the existence of 
entropy current with vanishing entropy production to determine the transport coefficients. 
We found all coefficients are uniquely determined up to 5 arbitrary functions, with the 
explicit expressions in (52) and (53). We also found that the stress tensor does not receive 
correction from terms proportional to the magnetic field. 

It will be interesting to extend the present analysis to chiral superfluid [1], which is 
a SU{Nf) X SU{Nf) SU{Nf) superfluid. Rich phenomena associated with anomaly 
effect on low energy finite density QCD have been discovered over the past few years, see 
e.g. [37, 38, 39]. It is certainly desirable to formulate a superfluid hydrodynamical theory, 
which will be useful for the exploration for the dynamical aspects of low energy QCD. 
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A Derivations of d^O'^ 

In this appendix, we want to express the divergence dfjiO^{0 = co, B, Q, K, U, V, L) in terms 
of our basis set: pseudovector- vector and pseudoscalar-scalar, by using the ideal hydrody- 
namical equations. 

In 3 + 1 dimensions, it is convenient to use vector notations. We will work in the local 
rest frame(LRF) of the normal component, where = (1,0). Note the condition = — 1 
gives dfj_u^ = 0, but d^u ^ 0. We start by solving the following ideal hydrodynamical 
equations: 



= 

Expressed in the vector notation, (54) takes the following form: 



(54) 



se + s = o 

wii + Vp + pVX -iC = nE 

c + f (e + V • e) + dtfe + v/2 . i= 

V// = - c vtr 

i+ = E 

Vxi= -B 

where the first equation is the conservation of entropy, which follows from a combina- 
tion of energy conservation and current conservation. The second and the third equations 
arc momentum and current conservation respectively. The fourth equation is from the 
"Josephson" equation. The last two equations arc vector form of the last equation of (54). 
u is readily solved as u = ^{nE — Vp — f^VX + ^Q, which will be repeatedly used below. 

In LRF, the gradient terms O'^ has the following simple expression: 



cJ = -V X «, B = -V X ^ 

n = ^ixit, K = -ixE 



(56) 
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where the quantities with over left arrows ^ V^T are to be contracted. The diver- 
gences of O'* are worked out in the following: 



2 

= ^{nE-u- ajf^di^p - /^w'^^^X + ^ ■ ooQ (57) 
di^B'^ = ^e'^'^^^a^u.F^^ = u-B-2E-u 

= ^{nE-B- Bi'd^p - fBi'd^X + ^-BQ-2E-uj (58) 

= -e^-co + fii-cj-^il-B- ujVfi -nd-u + ^e'^'^^jdtdiUk 

= -9^-Lo-f^d-uj + uj- Vvtr -^d-B-E uj + ^e'^'^Cjutui 

= -9^ -co- ^i-B^'daP - fB^duX + nE-B + ^-BC)- ^(nE ■ oj + ^ ■ ojC 
2w w 

- uji'd^p - foo'^df.X) + io'^eOf^u, + ^e'^'^'^u^^aUPdpd^u^ (59) 
d^K^ = ^d^P^'^e^^^peP''^ + \p'"'e,.apd^eF''^ - d^^fiB^ - fid^B^ 

= -e^ ■ B + ^^ii ■ B - ii ■ {i X E) + ^E ■ B - B -V n - lid ■ B + ^e'^'^^jdtdiAk 
= -e^-B + 2iiE-u + -{-Kf'd^p - fRf'd^X) + Ie-B- B^dali 
+ ^e^'^'^^u^^an'dpd^Af, (60) 
d^U'' = ^,f^--l^d^u,^ad^p+^e^'''''^u,d^^ad^p 

= ^u- {ix Vp) - iJ^Vp -^-up-^B -Vp 

= ^{nK'^d^p - 2fU^'d^X) - fiio^'d^p - e • uju'^d.p - ^B^'d^p (61) 
= ^{nK^di_,fi + 2Uf'd^ii - 2fV^d^X) - fiuj^d^fi - e • uu'^di.ii - \B>'d^ix 

= \(^ijk^Aiidkui + ]^u-{ix ^vtr) - iia\v\t -^-lu^-u-^b ■ Vvtr 

-^{nE-i- l^d^p - SH^d^X + /2() _ \B^^^d^u,. (62) 
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In (62) the first term is expressed as follows: 



eijk^jOi^iOkUi = ■ + s • I p 1- / -I I 

-U-Be + co- (VX + /x(V/x + 1"vtr)) + \a^J^B'' (63) 

B Relations among pseudovector vector 

The redundant pseudovector-vector are related to those chosen in (41) and pseudoscalar-scalar 
as follows: 

Q 

-l-{nE-uj- u^'d^p - fto^'d^X) (64) 

a^J^'K'' = K''C^^,u^ - E ■ ■ u + I'^E ■ u (65) 

Tf^J^U" = ~ (2L'*5^p - e • ojl'^df.p + l^u^'d/.p) (66) 

(T^^J^V = - ^ {2L^'d^fi - e • ojl^'d^fi + fu^'d^ix) (67) 

a^./'^L'^ = • ua^jn'' - ^^a^J^^u'' (68) 

^^^ed^u, = -l{^K^^-d^u, - L^d^p - fL^d^X) (69) 

U'^eO^^u, = -Li'd^p (70) 

Vi'C'd^.u, = -L^d^fi (71) 

L^Cd^u, = (72) 

^^d^P = --{^m^d^.u, - fU^d^X) (73) 

U^d^p = (74) 

V^'d^p = -V'd^ii (75) 

V^'^^^l = (76) 

n^'df.X = - ^ {nK^'^^X + 2U^'d^X) (77) 

^''d^^i = -^{nK^'d^ii + 2U^d^f, - 2fV^'d^X) (78) 
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F^,e^^ = '^i.J^B'' + ^-B-- Bf^C'd^u, (79) 

F^^,CB^ = (80) 

= -i-(n£;^ - d^p - f%X + ^^Cm ■ Be - i'B'^) (81) 

F^.CK^' = E-ii-B-fE-B (82) 

F^.eU'' = —I'^d^P^ ■ B + ^l^B^^d^p (83) 

F^uCV" = - Vs^/.^ • B + h^B^'d^fi (84) 

F^^eL'' = —I'F^^e^^^ - • Ba^Ji^r + ^fa^J^^B^ (85) 

E-^ = -^ {K'^d^p + fRf^d^X) (86) 

E-K = (87) 

S • [/ = ^if^a^p (88) 

£; . F = IjC'^a^// (89) 

E-L = ^Rf^ed^^u,. (90) 



Note that the scalar l^d^fi in (66) is not hsted as a independent scalar. By ideal hydrody- 
namical equations, it can be expressed as: 

i^di^fi = u^^d^x - fiu^d^fi + E-i—~^c + -{i^d^p + fi>^d^x) - a^jn^ - '^e. 

www f-f- 3 

(91) 
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